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Abstract 

The measurements of the 77 and bb partial decay widths of an intermediate- 
mass Standard Model Higgs boson are among the most important goals of a 

future photon linear collider. While in an initially polarized J z = state the 

2 

background process 77 — > bb is suppressed by — radiative corrections remove 
this suppression and are known to be very sizable at the one-loop level. In 
particular a new type of purely hard double logarithmic (DL) correction can 
even make the cross section negative at this order in perturbation theory. The 
second order term of the series is also known and enters with a positive sign at 
the two loop level. From a theoretical as well as a practical point of view it is 
clearly desirable to resum this series and to know its high energy behavior. In 
this paper, we derive the series of these novel "non-Sudakov" logarithms to all 
orders and calculate the high energy limit analytically. We also give explicit three 
loop corrections in the DL-approximation as a check of our results and show that 
the three loop structure also reveals the higher order behavior of mixed hard and 
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soft double logarithms. We are thus able to resum the virtual DL-corrections to 
77 {Jz = 0) — ► bb to all orders. 



1 Introduction 



Among the many important physics processes which can be studied at a future high- 
energy photon linear collider (PLC) [[I], [| , the production of Higgs bosons, 77 — > H, 
stands out. In particular, identification of an intermediate- mass Standard Model Higgs 
boson via its (dominant) H — > bb decay channels seems feasible and would provide an 
important measurement of the 77 and bb partial decay widths. 

The dominant (continuum) background to 77 —>■ H — > bb comes from the Standard 
Model processes 77 — > cc, bb. Various techniques can be used to suppress these. The 
most effective appears to be to polarize the initial-state photons: Higgs bosons are 
only produced from an initial state with J z = 0, whereas the leading-order qq (q = c, b) 
backgrounds are dominantly produced from the J z = ±2 initial state. More precisely, 
production from the J z = state is suppressed by m 2 q / s j|, |5|, |j. Hence in the region 
of the Higgs resonance o qq \J z = 0)/er 9f ~ m^/Mfj and the cc and bb backgrounds are 
heavily suppressed. 

Unfortunately, the 0{m 2 ) suppression of the J z = background does not sur- 
vive beyond leading order in QCD perturbation theory. In particular, there is no 
corresponding mass suppression of the 0(a s ) 77 — > qqg J z = amplitude squared. 
Although this gives rise notionally to three-jet configurations, in certain circumstances 
(for example, collinear 7 — > qq splitting followed by 7g — > qg Compton scattering) 
these can mimic the two-jet Higgs — > bb signal. Of course a combination of 6-tagging 
of both jets and three-jet veto can in principle eliminate the bulk of this background, 
and one is left considering the perturbative corrections from virtual and soft real gluon 
emission. These have the same 0(m 2 /s) suppression as the leading-order contribu- 
tionF] However the perturbation series now contains large double logarithmic correc- 
tions ~ a" \og 2n (s/m 2 ) at each order n, which potentially can give sizable corrections. 
Apart from their phenomenological importance, these logarithms are of theoretical in- 
terest. They do not have an obvious 'Sudakov' exponentiating behavior, arising from 
the fact that they originate in non-trivial regions of soft gluon phase space. Further- 
more, for any particular virtual gluon diagram they are mixed up with the genuine 
infra-red logarithms ~ a"log n (m^/A 2 ) log n (m^/s), where A is a fictitious gluon mass 
introduced to control the soft divergences. 

In Ref. the double logarithmic contributions were studied in considerable detail 
through two loops, i.e. up to and including 0(a 2 \og 4 (s/m 2 )) contributions. Inspection 
of the one- and two- loop coefficients led to the claim that, for m q = m b and y/s ~ Mh, 
the perturbation series was reasonably well approximated by the lowest two orders. 
However from both a theoretical and practical pont of view it is clearly desirable to 
have more information about the behavior of the double-logarithmic (DL) terms to all 
orders, in particular to see if resummation is possible. This is the goal of the present 

lT Jnlcss otherwise stated, from now on all statements about amplitudes and cross sections will refer 
to the J z = hclicity projection. 
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paper. As we shall see, a careful study of all the relevant two- and three-loop diagrams 
does reveal the pattern of the DL terms, and allows the all-orders result to be derived 
and resummed. 

In this paper we shall restrict ourselves to the theoretical study of the origin and 
calculation of the DL terms. In a future paper we will investigate their numerical effect 
on realistic PLC cross sections. We begin, in section 2, by deriving the leading order 
result for the J z = cross section, identifying the origin of the 0(m 2 ) suppression. 
In section 3 we discuss the structure of the one-loop corrections, show how the DL 
contributions arise, and compare these with the known one loop result ||, ^, [7J. The 
main part of the paper is section 4, where we calculate all three-loop DL contributions 
and deduce the all-orders result. Section 5 contains a summary of our results. 



2 The Born Amplitude 

We begin with a rederivation of the underlying Born process. Recall that for an 
initial state with J z = the cross section vanishes in the massless limit. In other 
words, the background to the standard Higgs production process via a fermion loop 

and subsequent decay into a bb or cc final state is suppressed by a factor of — . Higher 
order radiative (Bremsstrahlung) corrections remove this suppression in principle 0. 
One therefore needs to keep at least terms of O (77*) in the calculation of the Born 
amplitudes shown in Fig. [I]. We are also only interested in the region of large scattering 
angle 9, which is the dominant event topology of the Higgs signal || 0]. 

Keeping only terms of O y^j , we find it convenient to separate out the mass 
dependence of the spinors according to || 

u"(p) = £^\q-ti (1) 
«"(p) = (<?>-Ml Z fe= ( 2 ) 

V 2 p- q 

where the spinors on the r.h.s. are massless and q is an arbitrary reference four- 
momentum with q 2 = 0. The separation of the mass terms then allows for the conve- 
nient usage of the polarization vector || 
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k 2 , X 2 "P2 > 1^2 k 2 , X 2 "P2 j 1^2 



Figure 1: The Born process in polarized 7 + 7 — > 9 + 9 with the notation of momenta 
used in the remainder of this work. Keeping only the mass term in the numerator of 
the fermion propagator gives the correct Born result only for 9 ~ 0°. Details are given 
in the text. 



where here we denote the arbitrary massless reference momentum by h. In the following 
we will also use simply m for an arbitrary fermion mass, as the results derived do not 
depend on any heavy quark approximation. It is understood that for any application, 
the relevant fermion masses for the background to Higgs production below the W + W~ 
threshold are the b- and c-quark masses. The cross section of the latter is enhanced 
by a relative factor of 16 from the different charges of the quarks. 

Using the above notation we find for the Born amplitude M. Bom = B% + B 2 with 
an initial J z = state: 

_ 8me 2 Q 2 5 Al ,A2^i,A t2 ^Ai, At i 
Vs(l — p 2 cos 2 9) 



with (3 — yl — e the electromagnetic coupling with magnitude Q and 8 the 
scattering angle. The Kronecker-5's allow for a non-zero contribution only in case 
all four helicities are the same. The result agrees with Refs. [0, § up to O (^) 
and displays the somewhat surprising feature of a bosonic angular dependence (as in 
Bhabha scattering for instance), with no #-terms in the numerator. 

This feature leads one to suspect that by keeping only the mass term of the virtual 
fermion propagator, one should obtain Eq. |5| very easily. A detailed analysis of this 
contribution, however, yields: 

^Bom = Am t QH ^T [<*1, "All**, *l)(l + I 2 COS 2 g)<W + 

s(l — p cos 9) L 
(k h Ai|* a , -Ai)(l - (3 2 cos 2 0)<5_a 1iW ] (6) 

While it does not agree with the exact result in Eq. |5|, it is also evident that for 
small angles and high energies both solutions do agree. It was found in Ref. |7| that this 
feature holds also through one and two loops if one is interested only in the leading 
logarithmic accuracy at large angles. The reason for this seemingly contradictory 
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behavior is that while the angle between the outgoing fermion (antifermion) and the 
beamline is large, the angle between the incoming photon and the emitted soft parton 
is small. Therefore, while it is wrong to keep only the m-term in the propagator of 
the virtual fermion for the Born amplitude, it actually does give the correct leading 
logarithmic structure of higher order contributions with a soft fermion line. 

For higher order corrections it is also convenient to use the following expression for 
the J z = M.Bom amplitude: 

8me 2 Q 2 v 

Mbotu ~ r-r?:^ \(h, h) S2 XV {k 2 , h) (pi, X\p 2 , -A) (7) 

s sin v 

In the following we describe how to obtain the double logarithms in the one loop 
approximation and introduce the techniques applied in the later stages of the paper. 



3 One Loop Form Factors 

At the one loop level, the only diagram that can contribute large double logarithms 
is the box diagram with a virtual gluon connecting the outgoing qq pair. While the 
general method of employing Sudakov parametrizations and integrating over the per- 
pendicular components of soft loop momenta is well known |§ and has been used 



extensively |TI| [XT], [l^, [13] , we still find it useful to briefly review some basic proper- 
ties. In section f| we will also remark on some subtleties one encounters with higher 
order virtual corrections. 

Large logarithms occur only when there are small and large scales connected in 
a Feynman diagram. While the particular Feynman diagram is gauge dependent, the 
leading logarithmic contribution is not. In our case, the box diagram is connected with 
vertex and self energy corrections to form a gauge invariant cross section. However, in 
any gauge only the four-point diagram can yield the desired pair of large logarithms. 
It is therefore independent of the particular choice of gauge. A familiar and analogous 
example of the same mechanism is the gauge invariance of the infra red \ term in the 
context of dimensional regularization for a massless box diagram. 

In Fig. |2| we list the four distinct topologies that allow a doubly logarithmic 
contribution from the box diagram. The "blob" always denotes a hard momentum 
flowing through the omitted propagator and in each case there is a corresponding soft 
momentum flowing through the propagator facing the hard blob. 

The feature which makes this process interesting and very different from the stan- 
dard Sudakov form factor is the fact that we have soft fermion contributions as well. 
"Soft" means in this context soft compared to the hard blob. 

The Sudakov decomposition of loop momenta should be chosen such that the ex- 
ternal four momenta joining the soft parton line are used as the base vectors. For 

2 In higher orders this is not always unambiguous, as we will discuss in section ^ and can lead to 
technical difficulties. 
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Figure 2: The one loop soft and hard topologies contributing to the doubly logarithmic 
form factors. These graphs are obtained from the one loop box diagram. We denote 
the soft topology as A, the first hard topology as B and the second (third) as C. 



the graphs of Fig. [^, however, we take: 



2 \ / 2 

m \ I m 



topology A : h = fa ipt —p 2 I + a x I p 2 —pi I + h± 



topology B : I = (3k\ + ak 2 + l± (9) 

(771 2 \ 

P2 k 2 J + l± (10) 

At large angles and in the doubly logarithmic (DL) approximation, we have 

2px ■ ki = 2p 2 ■ k 2 w 2pi • k 2 = 2p 2 ■ k x « 2p x ■ p 2 w 2k\ ■ k 2 = s (11) 

The phase space for the large logarithms is given by 

1 > \a\,\/3\ > |Zl/s|,> \m 2 /s\ (12) 

1 > > \h 2 Js\ > |A 2 /s| (13) 

where we already have multiple gluon insertions in mind. The four dimensional 
Minkowski measure is rewritten according to d 4 l = d 2 l±d 2 l\\, with 

d 2 l± = \l±\d\l ± \d<t> = ^dl 2 ± d<f) = ndl 2 ± (14) 

d 2 l\\ = \d(l ,l x )/d(a,f3)\dadf3^^dadp (15) 

The integrations over the transverse momenta of the soft particles is performed by 
taking half of the residues in the corresponding propagators: 

dH s dad[3d 2 l± -2 s iin^( a ^ ( , a \ 

-m -da dp H(sap — m ) (lb) 



P — m 2 + is 2 sa(3 + l± — vn? + is 2 

S ^'^^.i -i^-da^^^saA-X 2 ) (17) 



If - A 2 + is 2 saipi + hi - A 2 + ie 2 



The DL-contribution of a particular Feynman diagram is thus given by 

M k = M B orn ?k (18) 

where the are given by integrals over the remaining Sudakov parameters at the 
n-loop level: 

^ ( ol s \ u ^ f 1 [ l da>i dfii ^ 
topology A : T k = — }[ / / — fc 

\2lXj tJ{J0 JO OLi jJi 

„ „ ^ (OL s \ n f l f 1 dad0 n 1 =} f 1 f l daiddi^ 

topology B , C : T k = M / / —^11/ / ~ 

\2irJ Jo Jo a p fj[ Jo Jo ot{ Pi 

The factors Q k are specific to the given Feynman diagram and contain its color 
factors and the restrictions on the occurring Sudakov variables necessary to ensure that 
the matrix element has logarithmic behavior in each of them. The range of integration 
is restricted by demanding that there is at least one pole term on the opposite side of 
the contour of integration from the other propagator terms, since Cauchy's Theorem 
would otherwise give a vanishing contribution. 



(19) 
(20) 



As a simple example, we now apply the above techniques to the one loop diagrams 
of Fig. [| For the standard infrared form factor of topology A the Sudakov phase 
space is schematically given in Fig. § and leads to the following result: 



T 



A 



-c F 
-c F 



a s r 1 r 1 da\ d0_ 
2n Jo Jo ai (3 



^ (sa l l - A 2 ) (at - ^-0^ U - 



2tt 



m dPi f 1 da\ f 1 dPi f 1 da% f~r d0% f 1 da± 







ati 



d(3\ r 1 da\ 



-a 



m 



m 



log h log — lo£ 



01 



A 5 



-01 Oli 



0i 



s/3i 



[2V 



2n \2 ~ s s ' m" 

The result is identical to the usual DL-Sudakov form factor 101. For the hard form 



factors of topologies B and C we have analogously: 



F = Tn = T C 



f 1 [ 1( ^ d J-Q(sa(3- 
Jo Jo a v 



da d0 



2w Jo Jo a 



m 



-C F — log — 
Air s 



(22) 




Figure 3: The one loop soft Sudakov phase space of topology A. The 0-functions 
given in Eq. ^ are translated into integrals covering the shaded area in the figure for 
the Sudakov parameters a\ and (3i as defined in Eq. |]. At higher orders, the phase 
space becomes increasingly untrackable and even numerical integrations need many 
iterations in order to converge. 

Taken together, and noting that topology C has a multiplicity factor 2, we thus 
have rederived the one loop structure of Refs. || |6| and in the DL-approximation. 
The fact that the Born amplitude factors out is trivial for topology A. For topologies 
B and C the exact numerator factorization was derived in Ref. 0. In the former, 
it is connected with the aforementioned purely bosonic angular dependence of M. Bom 
of Eq. ||, as it contains a hard subprocess q + q — > q + q. In the latter, one needs 
to exploit the Dirac equation and gauge invariance in order to factor out the hard 
subprocess 7 + q — ► q + g and the Born amplitude. Details are given in . 

With these results, we can now proceed to investigate the higher order structure 
of the leading logarithmic corrections from insertions of multiple gluons in the basic 
contributions shown in Fig. |2|. The next section gives explicit results through three 
loops which enable us to arrive at an all orders resummation of both hard as well as 
soft logarithmic corrections. 
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4 All Orders Resummation 



In this section we derive the expressions for the leading logarithmic corrections to all 
orders for the J z = initial state 7 + 7 — > q + q process. While it is possible to 
rigorously derive these corrections to all orders for topologies A and B, for topology 
C only the new hard "non-Sudakov" logarithms can be treated analytically. For soft 
insertions we derive explicit three loop results and show that essentially the same 
factorization structure emerges as in B. This behavior is then extrapolated to all 
orders. 

At this point we note that we have recalculated the virtual two loop corrections 
presented in Ref. |7j and agree with their results. In the following, we therefore restrict 
ourselves to a discussion of the three loop level. All resummed formulae were checked 
explicitly through one-, two- and three-loops. 



4.1 The Sudakov Form Factor 

We begin with the Sudakov form factor of topology A. At three loops, the relevant 



Feynman diagrams are given in Fig. 12. The summation is analogous to the soft 



vertex corrections where the exponentiation to all orders for a color singlet final state 



was proved in Ref. |15| to leading order by employing an asymptotic Hamiltonian as 
well as the coherent state formalism. Explicit three loop calculations were performed 
in the leading logarithmic approximation in Ref. [M and related work can also be 



found in Refs. |L7| |18| . For all orders corrections to topology A we thus find the 
familiar result: 

r oft = exp {T A ) (23) 

The non-Abelian color factors of diagrams with an "Abelian" topology are canceled 
by the contributions of the purely non-Abelian topologies. In this context it is worth 
noting that for the three gluon vertex coupling, there is an effective factor of — | 
stemming from two contributions in numerator terms that are needed to obtain the 
required number of logarithms. This is only the case if one has a denominator of the 
tyP e a-i-y ■ For one gluon very soft, there is a term —2U ■ p and one U ■ p, U being the 
other (relatively) soft gluon momentum and p an external hard fermion momentum. 
The Lorentz structure of the three gluon vertex is such that only these two terms 
give DL-contributions (the third connects the same fermion line). The effective soft 
contribution is therefore — I ■ p, compared with a factor 2p u from the emission of a soft 
gluon off a fermion line. 



4.2 Resummation of Novel Hard Logarithms 

The interesting feature of the higher order corrections to the one-loop diagrams of Fig. 
^| is that there are novel purely hard "non-Sudakov" logarithms due to the presence of 
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Asymptotic Hard Form Factor 




Figure 4: A comparison of the exact summed result of Eq. |4] with the asymptotic form 
given in Eq. and the first four terms in the expansion. The value of the coupling 
was kept fixed at a s — 0.1 and x = log At a fixed order in perturbation theory the 
high energy limit alternates. Only the exact and asymptotic form factors can be seen 
to give the correct limit at high energies. 



the "soft" fermion line in topologies B and C. In the case where a gluon insertion does 
not end on an external (hard) fermion line, the soft fermion line serves essentially as 
a cutoff for gluons that are soft compared to the hard blob. The diagrams are listed 



in Figs. [13] and [14] and denoted by hi, as they contain no soft Sudakov contributions. 
Besides a different color factor, both contributions are identical and lead to the same 
structure as in the soft case, except for the last two integrations. They have the 
additional effect of eliminating (a>i — ^yPi) terms, which are very important in the 
soft case (see Fig. The reason is simply that they are automatically fulfilled since 
the soft fermion line induces a — cutoff for the soft gluon Sudakov variables. The 
integrations over the gluonic variables can thus be performed easily and summing over 
i gluon insertions we have: 

r~ „ £ (|i)- + 7 1 a |/ 1 i ^(zip iog . /3log . a 

^(a s y +1 (-If /- 1 dp, i f m 2 V +1 
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Relative Contribution 




-0.3 

6 8 10 12 14 

x 

Figure 5: The relative contribution of the one (1L) and two (2L) loop corrections to the 
exact summed result of Eq. The value of the coupling was kept fixed at a s — 0.1 
and x = log For the b-quark mass, v^i = 100 GeV corresponds roughly to x = 6.2, 
y/s = 1 TeV to x = 10.8. In case of the c-quark, y/s = 100 GeV corresponds roughly 
to x = 8.4, y/s = 1 TeV to x = 13. The higher than second order corrections are thus 
important whenever the experimental accuracy reaches the percentile level. 

Using partial integration i-times (differentiating log l+1 integrating lo ^ - ) , we 
find: 
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(24) 
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In order to compare this result with the known one- and two-loop leading logarith- 
mic results we list them below explicitly and find agreement with Refs. || [| [7J: 

hard l (° i s\,- 2 m 2 1 fa s \ 2 4 m 2 1 f a s \ 3 6 m 2 



- ' • -(^) log 4 — -^y^l 1"K 




(25) 

The asymptotic behavior of the hypergeometric function 2-P2 for s — ► 00 can be 
obtained by using an integral representation and taking the leading pole term. In 
our case we find 

■hard a Sl ,m 2 /3\ r dt T 2 (l + t)T(-t) 



jrnara _ i \ z p 

47T 



'a, , n m 2 




ioo 27rir(2 + t)r(| + t 




- log — log 2 — (26) 
\n s 

The last line follows from the residue of the integral representation at t = — 1, 
closing the contour in the negative t-plane. All higher negative pole contributions are 
suppressed by at least a factor of — ^r- With the ^= suppression contained in the 

OS s 

Born amplitude @, the high energy limit is well behaved. 



It is also possible to reconstruct the asymptotic series based on the above integral 
representation. The higher order poles are simple ones and can thus be evaluated 
easily. We find 



t 



hard ^ _^ log2 ^! r ^ d^ T 2 (l + t)T(-t) /«. Tg\ 

4tt 6 S V2yi^oo2«r(2 + t)r(f + t) v 87r s / 

= -7 - log log 2 <*) + £ 2 " r (i) r (-^) (27) 

Although it possesses the ominous factorial growth of the coefficients in the series 
for the higher order residue contributions, the equivalent original expansion in terms 
of the confluent hypergeometric function in Eq. [24] means that it is Borel summable 
and that the high energy limit in Eq. ^ does indeed give the correct behavior for 
s — > 00. Fig. f| contains a comparison of the exact hard form factor of Eq. |24| and its 
asymptotic logarithmic divergence given by the first term in Eq. |27|. They can be seen 
to agree in the high energy limit. In addition the figure demonstrates that at a fixed 
order in perturbation theory, the result alternates between ±00 for high energies. 
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Asymptotic Residues 




20 25 30 35 40 45 



Figure 6: A ratio of the first four residues of the integral representation of Eq. |27] with 
the exact result of Eq. El. The value of the coupling was kept fixed at a s = 0.1 and 



x 



log 



At lower energies the 



T divergences of the higher order residues is 



clearly visible. At high energies, however, they can be seen to give the correct limit. 



The numerical effect of the higher order terms is shown in Fig. [| where the all 
orders hard form factor is compared to the one- and two-loop corrections at realistic 
future collider energies. As explained in the figure, only if the experimental accuracy 
reaches the percentile level are higher than second order contributions important. The 
relative size of the all orders correction is larger for m = m c than for the b-quark mass. 

A more detailed comparison of the behavior of the higher order residues is given 
in Fig. |6|. The first four pole contributions are plotted against the exact hard form 
factor and agree with Eq. |24] in the asymptotic regime. At lower energies the -, — 1 \ n 

behavior of the first three simple residues of Eq. |27| is clearly visible. 

In order to check this derivation explicitly at the three loop level, we have calculated 
the corresponding diagrams denoted by h in Figs. [13] and |14| and find the following 
results (modulo color factors and chosing for the Sudakov decomposition {^1,^2} for 
topology B and j& 2 , (j>2 — V'W)} f° r topology C): 

& hl ~ 6 (saP - m 2 ) 9 (a 2 - a x ) (a x - a) 9 (fa - A) 6 (A - P) (28) 
Q h2 ~ (sa/3 - m 2 ) { ai - a 2 ) 6 (a 2 - a) (/3 2 - /3 X ) Q (J3i - (3) (29) 
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Hard Form Factor 




Figure 7: A comparison of the third term of the series expansion in Eq. [24] (line) and 
the explicit three loop results given in Eqs. and (triangles) . For the three loop 
data we used 10 6 evaluations for each of the 100 iterations of the Monte Carlo integrator 
VEGAS [p0[| . nib = 4.5 GeV was used in both results. The absolute values are plotted 
and the agreement is below the percent level. Color factors and the coupling constant 
are omitted. 



The remaining terms h 3 and in each figure cancel the non-Abelian color factor 
stemming from the crossed box insertion h 2 , i.e. h 3 = /i 4 ~ -y. It can easily be 
seen from the explicit form in Eqs. |28] and |29| that the contributions of the gluon 

2 

insertions are cut off effectively by the fermion line Sudakov variables for which — is 
the minimum value attainable. 



Fig. [7| contains a comparison of the explicit three loop results given by hi + h,2 and 
the third term in the series [23|. The agreement is below the percentile level. 

For completeness, we give the form factors explicitly for the two topologies B and 
C below, including the coupling and their differing color structure: 

J%** = T 2 F 2 (1,1;2,~;~.F) (30) 
Jfrrd = jr aF2(ljl;2| 2 ; ^L^ (31) 

where JF is the one-loop hard form factor given in Eq. 
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Figure 8: A typical pair of three loop diagrams for topology B which vanishes in the 
DL-approximation when both contributions of a gluon radiated off a loop left of the 
hard blob couples to the outgoing fermion ( g% ) and the one coupling to the outgoing 
anti-fermion ( g 2 ) line are added. The origin of this cancellation, which proceeds to 
all orders, is the spinor line minus sign of the anti-fermion. Details are discussed in 
the text. 



4.3 Mixed Form Factors 

While the novel hard logarithms of the previous section are the most interesting new 
feature of the J z = 7 + 7 — > q + q process, technically the mixed hard and soft 
contributions are a considerable challenge. In fact for topology C we are only able to 
show numerically through three loops that the same structure emerges as for topology 
B. 



4.3.1 Topology B 

We therefore start with the diagrams given in Fig. O for topology B which contain a 
soft (s) contribution. The diagram S5 is evidently given by the product of the two-loop 
purely hard and the one-loop Sudakov form factors. Also for the sum of diagrams 
Si + ... + S4, it can easily be seen that we again have a factorization, this time the 
purely hard one-loop times the two-loop Sudakov form factor. 



The nice feature of this topology is that corrections connecting the original one- 
loop hard form factor with gluon insertions between the outgoing fermions vanish in 
pairs. The reason is simply a relative spinor line minus sign when a gluon starting 
from a fermion line of the original one-loop topology is connected with one of the 
final state fermions. One always has in this case an analogous correction which ends 
on the anti-fermion line, thus yielding a relative minus sign. Fig. || indicates the 
cancellation in pairs. In this particular case, for instance, we find (chosing the Sudakov 

parametrization (px ~ IT^i)} f° r 9i aric ^ {^i> (.P2 — V^i)} ^ or ^ 2 wr ^ n respect to 
the gluon going across the hard blob): 



S gi ~ -0 (saf3 - m 2 ) (sa^ - A 2 ) (sa 2 (3 2 - A 2 ) ( e*i - —ft 
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2 \ / 2 

m \ „ / „ m 



x e ipx - — aij e I & - — (« 2 - Pi) J e (« 2 - a) e (& - ai ) (32) 

6 S2 ~ 6 (aa/3 - m 2 ) 6 {sa x fa - A 2 ) 6 (sa 2/ 3 2 - A 2 ) \a x - ^-ftj 



2 \ / 2 



x I ft - — ai J e I /? 2 - — (ar 2 - ai) J 6 (a 2 - a) 6 (/3 2 - fa) (33) 

The g 2 contribution is minus that of g\ after the replacement ai < — > (3% has been 
made. The two diagrams therefore cancel each other. This feature holds to all orders 
and thus leads to a product of the hard form factor in Eq. |30] and the Sudakov form 
factor in Eq. In other words we have 

Tb = T 2 F 2 (1, 1; 2, |; X -T) [ exp (F A ) - 1] (34) 

where we have again used the notation introduced in section |3] for the one-loop form 
factors, hs refers to the mixed hard (h) and soft (s) nature of the form factor. 

4.3.2 Topology C 



For the diagrams given in Figs. |T4] and [15] there is no obvious way of generalizing 
a cancellation between diagrams to all orders. We therefore chose to calculate these 
corrections explicitly through three loops and compare with the "expected" analogous 
result found for topology B in Eq. |34|. The results are given in appendix 0. Fig. |] 
shows a comparison of the contribution from the second term of the hard form factor 
from Eq. |31] multiplied by the one-loop soft form factor of Eq. |21| with the three loop 
corrections given by s 13 + ...s 2 o i n Fig- HD- Only the diagrams Si 3 + Su + Si 5 actually 
contribute, the others are either zero (s% 9 and s 20 ) or cancel in the sum (si 6 + Si7 + Si8 = 
0) in the DL-approximation. The leading logarithmic behavior is shown to agree for 
both A — ► as well as for s — ■> oo, keeping m = fixed. The agreement is better 
than 1% and improves for large values of the logarithms. 



For the other mixed hard and soft contribution at the three loop level, namely the 
one-loop hard form factor in Eq. |22] multiplied by the second term in the Sudakov 
exponential of Eq. |23], many more diagrams contribute. The remaining soft diagrams 
of Figs. [14] and [15] are not only more numerous but also contain technical subtleties 
worth mentioning. 

In particular, diagrams containing a soft gluon radiating off the "soft" fermion 
line facing the hard blob need careful consideration. In those cases one has two (for 
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Figure 9: A comparison of the third term of the series expansion in Eq. |37| (line) with 
a double hard and single soft form factor (DhSs) and the explicit three loop results 
given in appendix |A.1.2| (triangles). The plot on the left keeps A 2 = 10~ 2 fixed, the 
one on the right has s = (100 GeV) 2 for each data point. For the three loop results we 
used 10 6 evaluations for each of the 100 iterations and m& = 4.5 GeV was used in both 
plots. The absolute values are plotted and the agreement is below the percent level. 
Color factors and the coupling constant are omitted. 



diagrams s n and s\2 of Fig. |TJ] even three) regions in the phase space contributing 
to the DL-corrections, corresponding to either side of the "soft" fermion line going on 
shell. 



For instance the two loop diagram in Fig. |10| possesses this ambiguity. It actually 
does not give a contribution in the DL-approximation due to an asymmetry of the 
diagram with respect to the external momenta {^2,^2}; while the large variables (see 
Eq. |TT| ) do not differ in the DL-approximation if we make the replacement k<i < — ► 
The three loop diagrams s±g and S20 in Fig. O vanishes for the same reason. The 



easiest way to see this on a technical level is therefore to parametrize l\ choosing 
{k 2 ,P2} for the Sudakov decomposition. Expanding around the left side of the "soft" 
fermion line in Fig. [Hi we find: 
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Figure 10: A two loop diagram for topology C which vanishes in the DL-approximation 
when both contributions of the soft fermion momentum I to the left and right of the 
gluon insertion are treated properly. As there is no "obvious" choice for the Sudakov 
parametrization of the gluon momentum, technical difficulties can arise. Details are 
discussed in the text. 



0i ~ ©(sa/3-m 2 ) [0 - «i) (/3 - /3i) («i/3 - a/3i - a/3 - «i/3i) 

+ (ai - /3i) S(J3- Pi) (a/3i - a x p - ap - a^fr)] (35) 

Shifting the fermion loop momentum, effectively I — > l+h, and expanding around 
the right side of the gluon insertion gives analogously: 

2 ~ —0 (sa(3 — m 2 ) [0 (a,\ — (3\) (a — a{) {aj3\ — a±(3 — a(3 — a\(5i) 

+ (Pi — «i) (a — Qi) (aiP — a>Pi — aP — QiPi)] (36) 

where we have omitted the common color factor. The second contribution is evidently 
(— ) the first after the replacements a < — > P and cti < — > Pi. For the soft contributions 
of the three loop diagrams S5, Sq, sq, sio, sn and su where the antisymmetry is removed 
by the additional gluon, we find it technically more convenient to choose {&2,pi} for 
the Sudakov decomposition for the genuine DL-corrections listed in appendix 0. A 
careful treatment should in any case give the same result, however. 

Fig. |TT] shows that the three loop contribution of the sum of diagrams Si + ... + S12 
is indeed equal to the product of the one-loop hard form factor in Eq. multiplied 
by the second term in the Sudakov exponential of Eq. The leading logarithmic 
behavior is shown to agree for both A — > as well as for s — ► 00, keeping m = m;, 
fixed. The relative difference actually decreases with increasing values of the large 
logarithms. 

It is therefore more than suggestive that the structure which prevails for the one-, 
two- and three-loop DL-contributions exponentiates to all orders analogously to that 
of topology B in Eq. 0, namely that 

T h c s = T a F 2 (l, 1; 2, ?; [ exp (f A ) - 1] (37) 
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DsSh Form Factor 




Figure 11: A comparison of the third term of the series expansion in Eq. [37| (line) with 
a double soft and single hard form factor (DsSh) and the explicit three loop results 
given in appendix |A.1.1| (triangles). The plot on the left keeps A 2 = 10~ 2 fixed, the 
one on the right has s = (100 GeV) 2 for each data point. For the three loop results we 
used 10 6 evaluations for each of the 100 iterations and m& = 4.5 GeV was used in both 
plots. The absolute values are plotted and the agreement is below the percent level. 
Color factors and the coupling constant are omitted. 

5 Summary 

In this paper we have derived a novel series of hard double logarithms in polarized 
11 {J z — 0) collisions to all orders. The resulting series is given by a confluent hyper- 
geometric function 2 -^2 which possesses a log log 2 ^j-^j high energy limit. Taking 
into account the ^= suppression contained in the Born amplitude, the new corrections 
are well behaved as s — ► oo. The series expansion agrees with the known one and 
two loop results as well as our explicit three loop calculation. 

It was furthermore shown that at three loops all additional doubly logarithmic 
corrections are described by the exponentiation of Sudakov logarithms at each order 
in the expansion of the new hard form factor. This behavior is already present at one 
and two loops and can thus safely be extrapolated to all orders. 
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In summary we list the complete virtual DL-form factor contribution to the process 
7 + 7 (Jz = 0) — > q + q: 

M DL = M B om {exp {F A ) + T 2 F 2 (1, 1; 2, | + 2 JF 2 F 2 (1, 1; 2, | ^-.F) 
+JF 2 F 2 (1,1;2,|^) [exp (^)-l] 

+2 F 2 F 2 (1, 1; 2, | ^-F) [ exp {T A ) - 1]} (38) 



Bom 



T 2 F 2 (1, 1; 2, | ~F) + 2 F 2 F 2 (1, 1; 2, ~; ^^)| exp {T A ) 



where T A and denote the soft and hard one-loop form factors of Eqs. and ^2 
respectively. It should be pointed out that the contributions containing the soft form 
factor T A cancel in a physical cross section with the infrared divergent contributions 
from the emission of real soft gluons. For realistic future collider applications Fig. 
[5] demonstrated that the higher than second order contributions are important if an 
accuracy in the percentile regime could be reached. A more thorough study of the 
phenomenological implications of this work including realistic experimental cuts is in 
progress . 
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A Three Loop Results in the DL- Approximation 



A.l Hard and Soft Corrections 

In this appendix we give explicit three loop results for the terms of Eq. 



of 



14| with the corresponding 
T3] the left and right 



topology C. The particular contributions are listed in Fig. 
notation. For the three loop contribution of topology B in Fig. 
side of the hard blob factorize for the graphs s\, S5. On a technical level they are 
therefore effectively two loop corrections. See Ref. for explicit results. 

In the following we distinguish between corrections with one and those with two 
purely hard form factors. We also omit color factors at this point, as the exponentiation 
of the soft gluons ensures that all "higher" order Ca contributions are canceled by 
purely non-Abelian corrections [0, [15], IT]. The correct color factors of the sum of the 



graphs given in Figs. 113, O and [16] are contained in the results of section 
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A. 1.1 Double Soft Single Hard Contributions 



The following results for the 0& correspond to the respective graphs of Fig. |I4| with two 
potentially soft gluons. The Sudakov parametrization for the three soft momenta (as 
indicated in the figure for si) was chosen to be jfc 2 , (l>2 ~ ^7-^2)} for the soft fermion 
line I in all graphs, j [pi — ^p-i) , (p2 — ^Pi) } fo r h in Si, s 2 and for l 2 in si, s 6 
and |/c 2 , (pi - ^-^2)} for h in s 3 , ...Si 2 and for Z 2 in s 7 , ...s 12 . 

e (« 2 - ^a) e (a - ^a 2 J e (« 2 - ai ) e (a - a) e (/3 - a) (39) 

6 S2 ~ -6 - ^pj e (a 1 p 1 - ^pj e ( a 2 p 2 - ?pj 9 (an - ^a) 

9 (p t - ^a 2 ^j (a 2 - ^AJ 9 (a - ^-a 2 ^j (a 2 - ai) (A - A) 
©09 -A) (40) 

e S3 ~ -e («/3 - ^ e («iA - 7J © (a 2 A - 7) © («i - ^-a) 

9 (fa - "p-a^ © (a - ^ (« 2 + A)) © («2 - a) © (A - an) 6 (/? - A) (41) 

e S4 e (ai3 - ?pj e (^a - 7) © ("2/^2 - 7) © - 

(771^ \ / 777.2 \ 

A - — «i J © I ai - — (a 2 + A) J © (« 2 - a) (ai - A) © {P - A) (42) 

(777 2 \ 
«i A I 6 (a 2 — «i) 6 ((a + /3) A — aa 2 — a/5 — a 2 A) 

a 2 - — (A + AO ){©(«- A) © (P - A) -©(/?- A - A - « 2 )} (43) 
e S6 (ap - ^pj e ^a - 7) © (^A - 7) © (ft - 7^1) 

© (A — — "2 1 © (on - A) © ((a + A a 2 - "A - a/5 - a 2 A) 
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on- — (Pi + a 2 ) J {6 (a - a 2 ) (P - ft) - 6 (P - ft - ft - a 2 )} (44) 

e S7 e (a(3 - ?pj e - ^ e (ft - 

© (ft - ^«2^) © (« 2 - «i) (ft - ft) © («i - «) (45) 
S8 © (a/3 - ^ (aift - ?pj (a 2 ft - ^ (ft - ^a^ 

(ft - ^e*i) (a 2 - a) (ft - ft) (a x - a 2 ) (46) 

S9 ~ -0 (a/3 - ^pj (aift - ?pj (a 2 ft - ^ (ft - ^a 2 ^j 

Pi- — ("i + « 2 ) ) © ((a + P)a 2 - ap 2 -a/3- a 2 p 2 ) 
0(ai -a){0(a-a 2 ) -0(/3-ft -ai)} (47) 

510 © (a/3 - ^ (aift - ?pj (a 2 ft - ^ (ft - (P 2 - pi) 

(TV? \ 
ft — («i + a 2 ) J ((a + P) a 2 - ap 2 — aft — a 2 ft 2 ) 

Q(a 1 -a){e(a-a 2 )-e(P- p 2 -a 2 )} (48) 

511 © (aft - ?pj (aift " y ) © (a 2 ft - y ) © (ft - (ft - ft) 

ft - — (ai + a 2 )j 

{0 (/3 - ft - «i - ft - a 2 ) ((a + ft a 2 - aft - a/3 - a 2 ft) 
((a + ft ai — (a + ft a 2 — aft — aift — a 2 ft — aift) — 
(a — a 2 ) (/3 — ft — ai) ((a + ft ai — aft — aP — aift) 
((a + ft a 2 - aft - a/3 - a 2 ft) + 

(a - ai - a 2 ) ((a + ft) a 2 - (a + ft ai - aft - aift - a 2 ft - a 2 ft) 
((a + ft ai -aft -a/5 -aift)} (49) 

512 (a/3 - ?pj (^ft - y) 6 (a 2 ft - ^ (ft - ^a 2 ) (ft - ft) 

( 777^ \ 
ft - — (ai + a 2 )j 

{0 (ft — ft — ai — ft — a 2 ) ((a + ft ai — aft — aft — aift) 
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6 ((a + (3) a 2 - (a + (3) a% - aft - aift - a 2 ft - a 2 ft) - 
9 (a — «i) 9 (/3 — ft — 0:2) 9 ((a + (3) a.\ — aft — a(3 — aift) 
9 ((a + (3) a 2 - a(3 2 - a/3 - a 2 ft) + 

9 (a — ai — a 2 ) 9 ((a + ft a x — (a + ft a 2 — a(3i — aift — a 2 (3\ — aift) 
9 ((a + (3) a 2 - a(3 2 - a(3 - a 2 ft)} (50) 

A. 1.2 Double Hard Single Soft Contributions 

Below we give the results for the 9& corresponding to the respective graphs of Fig. 



14| with one potentially soft gluon. The Sudakov parametrization for the three soft 
momenta (as indicated in the figure for s{) was {/c 2 , (p 2 — ^^2)} f° r the soft fermion 

line I and for li in all graphs and for l 2 in Si 6 , s 20 , {{pi ~ T"^) > (p% ~ V'P 1 )} ^ or 
l 2 in S13 and {fc 2 , (pi - ^2)} for l 2 in s M , s i5 . 

9 513 9 L(3 - ^-j 9 (a 2 ft - ^ 9 (a 2 - ^ft^) 9 (ft - 

9 (/3 - ft) 9 (a x - a) 9 (ft - ft (51) 

9 514 9 (a/3 - 9 (a 2 ft - ^ 9 (ft - ^a 2 ) 

9 (ai - a 2 ) 9 (a 2 - a) 9 (A - ft (52) 

9 515 9 L(3 - 9 (a 2 ft - ^ 9 (ft - ^a 2 ^) 

9 (a 2 - ai) 9 (a x - a) 9 (A - ft (53) 

9 516 9 (a(3 - ^-j 9 (a 2 ft - ^ 9 (ft - ^a 2 ) 

9 (ai - a) 9 (a 2 - a) 9 (A - ft 9 (/3 - (3 2 ) (54) 

9 517 ~ e(a/3-^e(a 2 /3 2 -^je(3 2 -^a^ 

9 (ai - a 2 ) 9 (a 2 - a) 9 (ft - ft 9 (/? - (3 2 ) (55) 

9 518 ~ e^ af 3-^&(a 2 f3 2 -^&(^ 2 -^a 2 ^ 

9 (a 2 - ai) 9 {a x - a) 9 (ft - ft 9 (/? - ft) (56) 

9 519 ~ (57) 

9 520 ~ (58) 
The last two diagrams vanish in the DL-approximation for the same reason as the 



two loop contribution given in Fig. |10[ The diagrams s\§, Si§ cancel each other after 
writing 
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a 



e 



«2 



a 



) = («i — tt 2 ) © («2 — Oi) + (Q!2 — cui) («i 



a: 



(59) 



Strictly speaking, the above cancellation takes place only for the C\Ca terms of 
amplitude s\%. The CpC\ term, absent in graphs Si 6 and s\t of Fig. 0, is canceled 
by the analogous contribution of graphs S21 and S22 of Fig. 115. 



In the same way the above mentioned cancellation of "higher order" Ca terms takes 

The remaining diagrams of topology C 



place with all non-Abelian graphs of Fig. [15 
in Fig 



16 all cancel among themselves in groups of two or three. 
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Figure 12: The three loop Feynman diagrams contributing to the purely soft form 
factor of the topology A. The sum of all terms exponentiates as usual with the non- 
Abelian diagrams canceling the Ca parts of the crossed "Abelian" graphs. 
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Figure 13: The three loop Feynman diagrams contributing to the purely hard (h) and 
soft (s) form factors of the topology B. The sum of all soft terms exponentiates from 
each order of the hard process. The remaining diagrams cancel pairwise as described 
in section [|. 
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Figure 14: The three loop Feynman diagrams contributing to the hard (h) and soft 
(s) form factors of the topology C. The sum of all soft terms exponentiates from each 
order of the hard process. 
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Figure 16: The three loop Feynman diagrams of the topology C which cancel either 
pairwise or in groups of three in the leading logarithmic approximation. 
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